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A  Test  of  the  Climatic  Normality  of  Synthetic 
Profiles  of  Atmospheric  Liquid  Water 
Content  Over  the  USSR  in  1973 


1.  INTRODUCTION 

The  stuff  that  clouds  are  made  of  can  be  a  decisive  factor  in  processes  of 
military  importance;  for  instance,  the  erosion  of  exposed  surfaces  on  high-speed 
vehicles  and  the  attenuation  of  electromagnetic  radiation  along  a  path  through  the 
atmosphere.  The  mass  concentration  of  hydrometeors  is  a  property  of  frequent 
concern,  but  size  distribution  and  other  parameters  also  play  a  significant  role  at 
times.  The  term  "hydrometeor"  encompasses  water  particles  of  all  sizes,  the 
large  ones  called  precipitation,  as  well  as  the  smaller  cloud  particles. 

Following  a  frequent  convention,  the  mass  concentration  of  hydrometeors, 
including  snow  and  ice  in  all  of  its  forms,  will  be  denoted  liquid  water  content 
(LWC)  throughout  this  report.  The  imprecise  qualifier  "liquid"  serves  as  a 
reminder  that  what  is  not  included  in  LWC  is  water  in  the  gaseous  phase,  which 
is  almost  invariably  the  more  massive. 

LWC  and  the  other  microphysical  properties  of  hydrometeors  are  not 
observed  routinely  by  the  weather  services  of  the  world;  consequently,  there  are 
no  archives  on  which  to  base  climatologies  of  these  parameters.  Instead,  short 
of  instituting  a  special  observational  program,  one  can  generate  such  a  clima¬ 
tology  only  indirectly,  through  the  use  of  correlations  established  between  the 
hydrometeoric  parameter  and  weather  parameters  that  are  observed  routinely 

(Received  for  publication  10  July  1979) 
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and  archived.  The  latter  was  the  approach  of  the  Environmental  Definition  Pro¬ 
gram  (EDP)  conducted  during  the  years  1973-1977  by  the  Air  Force  Geophysics 
Laboratory  (AFGL).  One  of  the  more  ambitious  endeavors  of  this  program  was 
the  "ll -station  study"  based  on  AFGL's  second-generation  EDP  "model"  com¬ 
monly  called  AFGL-2.  * 

AFGL-2  is  a  tedious  procedure  which  starts  with  a  manually  analyzed  vertical 
time-section  (altitude  vs  time)  incorporating  all  meteorological  data  that  are  sys¬ 
tematically  available  for  the  particular  location.  Then,  wherever  clouds  or  precipi¬ 
tation  are  reckoned  to  have  existed,  values  of  LWC  are  assigned  according  to 
typical  values  for  the  class  of  hydrometeor,  temperature,  and  other  circum¬ 
stances.  The  11 -station  study  applied  this  procedure  for  the  12 -month  period 
starting  1  February  1973  to  each  of  11  meteorological  stations  in  the  USSR.  These 
stations  that  were  selected  for  the  supposed  diversity  of  their  LWC  climatology, 
as  well  as  for  the  availability  of  suitable  meteorological  records,  are  identified 
in  Table  1,  and  their  locations  are  shown  in  Figure  1.  For  each  of  these  stations, 
and  for  every  3  hours  of  the  entire  year,  a  vertical  profile  of  LWC  was  extracted 
from  the  time-section.  Thus  the  ultimate  product  of  that  study  was  some  32,  000 
profiles  of  LWC. 

A  natural  question  to  ask  is.  Was  this  year,  on  which  such  analytical  effort 
was  lavished,  anything  like  a  typical  year?  An  early  attack  on  this  question  was 
made  by  Lenhard1  who  compared  1973  with  a  10-year  period  in  terms  of  precipi¬ 
tation  and  found  1973  somewhat  drier.  A  more  direct  attack  is  the  subject  of  this 
report. 

2.  THE  TEST 

Because  of  the  magnitude  of  effort  that  would  have  been  needed,  AFGL-2  was 
not  used  in  this  LWC-comparison  of  1973  with  a  longer  span  of  time;  AFGL-1  was 
employed  instead.  AFGL-1  is  another  model  for  estimating  LWC  from  conven¬ 
tional  meteorological  data.  It  is  a  decision  tree  that  uses  a  radiosounding  and  the 
surface  observation  as  input.  An  automated  version  was  used  for  this  study. 

Unlike  AFGL-2  which  treats  both  cloud  and  precipitation  LWC,  AFGL-1  prescribes 


T!  Lenhard,  Robert  W.  (1974)  The  Use  of  1973  as  a  Sample  of  Erosion  Clima¬ 
tology  over  USSR.  INAP  No.  98,  Internal  Report  of  the  Design  Climatology 
Branch,  Aeronomy  Laboratory,  AFCRL,  56  pp. 

*  Prior  to  1976,  what  is  now  AFGL  was  part  of  the  Air  Force  Cambridge 
Research  Laboratories  (AFCRL),  which  no  longer  exists.  Previously, 
AFGL-2  was  known  as  AFCRL-2. 
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Table  1.  The  Eleven  USSR  Stations  of  the  Environmental  Definition  Program 


Station  Number 

Name 

Latitude 

Longitude 

221130 

Murmansk 

68°  58’  N 

33°03'  E 

260630 

Leningrad 

59° 58'  N 

30° 18'  E 

276120 

Moscow 

55° 58'  N 

37° 25'  E 

333450 

Kiev 

50°  24'  N 

30°  27 '  E 

339460 

Simferopol 

45°01 '  N 

33° 59'  E 

282250 

Perm 

58°  01 '  N 

56° 18'  E 

352290 

Aktyubinsk 

50°  20'  N 

57° 13'  E 

361770 

Semipalatinsk 

50° 21'  N 

80°  15'  E 

384570 

Tashkent 

41°  16'  N 

69°  16'  E 

307580 

Chita 

52° 01'  N 

113°  19'  E 

315100 

Blagoveschensk 

50°  16'  N 

1 27°  30'  E 

the  water  content  of  only  the  larger  precipitation  particles.  By  definition,  all 
forms  of  ice  are  considered  to  be  precipitation  and  are,  therefore,  covered  by 
AFGL-1. 

Although  an  attempt  has  been  made  to  assess  the  absolute  accuracy  of 
2 

AFGL-1  and  AFGL-2  the  results  are  indecisive,  due  to  the  rather  few  cases 
available  in  which  there  was  an  independent  measurement  of  LWC.  The  accuracy 
of  AFGL-1  is  immaterial,  however,  to  present  purposes,  which  require  only  that 
any  bias  in  AFGL-1  be  the  same  for  1973  as  for  the  longer  interval.  At  present 
this  is  an  unverified  assumption. 

Using  AFGL-1  for  the  test  of  climatic  normality  was  facilitated  by  the  pre¬ 
existence  of  6-hourly  AFGL-1  profiles  of  LWC  for  the  11  USSR  stations  for  Febru¬ 
ary,  May,  July,  and  October  1973.  These  had  been  generated  for  an  intercom¬ 
parison  of  AFGL-1,  AFGL-2,  and  yet  another  LWC  model.  ^ 

The  longer  period  selected  for  comparison  with  1973  was  the  5-year  period 
from  1968-1972.  For  each  of  the  4  seasons,  randomly  selected  cases  from  the 
0000Z  and  the  1200Z  radiosoundings  comprised  the  5-year  sample.  The  seasons 
were  defined  as:  winter,  December-February;  spring,  March-May;  summer, 
June-August;  autumn,  September-November.  (The  reason  for  the  particular 

2.  Peirce,  Russell  M. ,  Lenhard,  Robert  W. ,  and  Weiss,  Bernard  F.  (1975) 

Comparison  Study  of  Models  Used  to  Prescribe  Hvdrometeor  Water  Content 
Values.  Part  1:  Preliminary  Results,  AFCRL-TR-7 5-0470,  20  pp. 

3.  Touart,  C.  N.  ,  and  Izumi,  Y.  (1979)  Comparison  Study  of  Models  Used  to 

Prescribe  Hvdrometeor  Water  Content  Values.  Part  II:  USSR  Data. 
AFGL-TR-79  (in  press). 


choice  of  months  in  1973,  or  of  "matching  seasons  for  the  5-year  sample,  is  lost 
in  the  unrecorded  past.  ) 

in  the  random  selection,  any  case  for  which  the  input  data  were  found  defective 
was  rejected,  and  the  process  continued  for  each  season  until  100  cases  had 
accumulated  for  which  the  total  sky  coverage  was  fi/8ths  or  more.  An  I.WC  pro¬ 
file  was  then  generated  for  each  of  these  cases.  Meanwhile,  the  cases  which  had 
been  selected,  but  for  which  the  sky  coverage  was  less  than  6/8ths,  were  arbi¬ 
trarily  assigned  an  LWC  of  0  at  all  altitudes.  This  black -or- white  treatment  of 
the  profile,  depending  on  the  sky'  coverage,  is  a  crude  device  for  factoring  into 
the  statistics  the  probability  that  a  randomly  positioned  trajectory  will  encounter 
the  particular  LWC  profile.  This  is  a  common  practice,  although  4/8ths  sky 
cover  is  the  more  frequent  threshold  value. 

The  individual  station-month  samples  for  1973  range  in  size  from  64  to  123 
profiles,  except  for  one  station-month  that  is  missing  altogether.  The  individual 
station-season  samples  for  the  5-year  period  range  in  size  from  115  to  254  pro¬ 
files.  The  grand  totals  are  4343  and  6372,  respectively,  for  1973  and  the  5-year 
period. 

The  comparison  is  made  in  terms  of  two  integrals  of  the  profile  of  LWC: 
integrated  LWC  (ILWC)  and  the  Environmental  Severity  Index  (ESI). 

ILWC  s S"  M  dz. 

0  1 

ESI  =/  °°  z„M  dz„ 

o  Z  Z 

where 

M  is  LWC  (g  m'3) 

Zj  is  height  (m) 

Zg  is  height  (km) 
so  that 

ILWC  is  in  kg  m~2 

ESI  is  in  g  km2  m  3 

(ESI  is  a  parameter  of  particular  significance  in  the  erosion  of  ballistic  reentry 
vehicles. ) 

To  facilitate  the  comparison,  cumulative  distributions  were  constructed  for 
each  season,  for  each  station,  and  for  both  integrals.  From  these,  certain  sta¬ 
tistics  were  evaluated:  average,  median,  90th  percentile,  and  95th  percentile. 

(The  90th  percentile  is  the  value  below  which  90%  of  the  sample  lies.  )  Because 
the  distribution  functions  tended  to  be  rough,  the  median  and  percentile  values 
were  not  derived  by  simple  interpolation  of  points  on  the  curve.  Instead,  inter¬ 
polation  was  made  via  a  quadratic  function  fitted  by  least-squares  to  neighboring 
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points  of  the  distribution.  The  number  of  points  used  was  decided  subjectively, 
case-by-case,  and  varied  from  3-10. 


3.  RESULTS 

The  basic  distributions  and  detailed  statistics  are  presented  in  the  Appendix. 

Figures  A-l  through  A-ll  compare  1973  and  the  5-year  sample  in  terms  of 
the  cumulative  frequency  distribution  of  ESI  by  season  for  each  of  the  11  stations. 
Figures  A-l 2  through  A-22  make  the  same  comparison  for  ILWC. 

Figures  A-23  through  A-25  are  bar  graphs  comparing  the  two  samples  by 
station  and  season,  in  terms  of  average  ESI  and  ILWC. 

Tables  A-l  through  A-4  make  a  numerical  comparison  by  station  and  by  sea¬ 
son  in  terms  of  average  ESI,  median  ESI,  90th  and  95th  percentiles  of  ESI.  A  is 
positive  when  73  exceeds  5-yr.  Tables  A-5  through  A-8  are  the  analogues  for 
ILWC. 

Table  2  compares  the  gross  averages  of  the  statistics  for  1973  and  5-yr.  By 
way  of  explanation,  the  first  entry  in  this  tabulation  was  derived  from  Table  A-l. 
All  43  values  of  A  were  averaged,  and  this  average  was  then  divided  by  the  average 
of  the  43  73-values.  The  result  as  a  percentage  is  the  entry  in  Table  2. 


Table  2.  Difference  between  1973  and  5-yr,  Averaged 
Over  All  Stations  and  All  Seasons,  Expressed  as  a 
Fraction  of  the  Corresponding  Average  for  1973 


Statistic 

Average  Difference 
1973-(5-yr) 

ESI:  AVERAGE 

2% 

MEDIAN 

-5 

90TH  PERCENTILE 

7 

95TH  PERCENTILE 

2 

ILWC:  AVERAGE 

-6 

MEDIAN 

5 

90TII  PERCENTILE 

-4 

95TH  PERCENTILE 

-7 
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1.  CONCLUSIONS 


A  mere  scan  of  the  figures  is  sufficient  to  sow  the  suspicion  that  there  is  no 
great  systematic  difference  between  the  two  samples  overall.  True,  there  are 
local  exceptions,  notably  Tashkent  where  1973  is  drier  throughout. 

The  suggestion  that  these  samples  are  cut  from  the  same  bolt  is  reinforced 
by  Table  2,  which  shows  that  the  differences,  when  averaged  over  all  stations  and 
all  seasons,  are  quite  small  for  all  4  statistics  and  for  both  integrals.  Also, 
these  differences  are  evenly  balanced  between  plus  and  minus. 

To  guard  against  the  possibility  that  this  averaging  may  smother  significant 
details,  a  test  was  made  of  the  plausibility  that  the  observed  local  differences 
could  have  arisen  purely  by  chance.  How  the  test  works  will  be  illustrated  on  the 
fifth  entry  of  Table  3.  This  entry  was  derived  from  Table  A-5  by  first  deleting  all 
As  which  are  less  than  10%  of  the  corresponding  73  value.  These  cases  are  arbi¬ 
trarily  judged  to  be  ambiguous.  Then  a  simple  count  of  -+'s  and  -'s  was  made  among 
the  residual  As.  The  result  in  this  instance  is  15  +'s  and  21  -'s.  Under  the  null 
hypothesis  that  the  1973  and  the  5-year  sample  come  from  the  same  parent  popula¬ 
tion,  the  +'s  and  -'s  are  equally  likely,  and  the  binomial  distribution  gives  8.  1%  as 
the  likelihood  of  a  15/21  split  in  36  Bernoulli  trails. 

The  same  procedure  was  applied  to  each  of  the  statistics  of  both  integrals, 
even  though  these  are  not  truly  independent.  The  results  are  summarized  in 
Table  3.  Except  in  two  of  the  eight  cases,  the  chance  probability  is  comfortably 
high,  giving  no  reason  to  suspect  any  fundamental  difference  between  the  two  sam¬ 
ples.  For  comparison,  a  50/ 50  split  of  +'s  and  -'s,  which  is  the  most  likely  out¬ 
come,  has  less  than  a  14%  chance  of  occurrence  in  the  circumstances  of  Table  3. 

In  view  of  all  of  the  foregoing,  the  final  conclusion  must  be  that  this  study  has 
uncovered  no  reason  to  doubt  that  1973  was  a  perfectly  representative  year  with 
respect  to  LWC  over  the  11  stations.  Of  course,  this  judgment  is  based  on  a 
comparison  with  a  small  5-year  sample,  but  whether  that  sample  is  itself  climati¬ 
cally  representative  will  be  left  as  an  exercise  for  the  reader. 
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Table  3.  Pure  Chance  Probability  of  the  Observed 
Combination  of  Differences  between  1973  and  5-yr 
for  Individual  Stations  and  Seasons 


Statistic 

Chance  Probability  of 
Observed  Difference 

ESI:  AVERAGE 

12.9% 

MEDIAN 

0.6 

90TH  PERCENTILE 

11.8 

95TH  PERCENTILE 

12.  6 

ILWC:  AVERAGE 

8.  1 

MEDIAN 

6.7 

90TH  PERCENTILE 

10.  8 

95TH  PERCENTILE 

3.4 
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Cumulative  Frequency  Distribution  of  ESI  for  Station  221130,  MURMANSK 


Figure  A2.  Cumulative  Frequency  Distribution  of  ESI  for  Station  260630,  LENINGRAD 


Figure  A5.  Cumulative  Frequency  Distribut 
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Figure  All.  Cumulative  Frequency  Distribution  of  ESI  for  Station  315100,  BLAGOVESCHENSK 
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Cumulative  Frequency  Distribution  of  ILWC  for  Station  221130,  MURMANSK 


leninbr 


ncy  Distribution  of  ILWC  for  Station  260630,  LENINGRAD 


INTEGRATED  LWC  C  KB/M*  1  INTEGRATED  LWC  t  KB/M*  3 

Figure  A14.  Cumulative  Frequency  Distribution  of  ILWC  for  Station  276120,  MOSCOW 


Figure  A15.  Cumulative  Frequency  Distribution  of  IL,WC  for  Station  333450,  KIEV 
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Cumulative  Frequency  Distribution  of  ILWC  for  Station  339460,  SIMFEROPOL 


WINTER 


roquoncy  Distribution  of  IIAVC  for  Station  2822.10,  PKKM 
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Cumulative  Frequency  Distribution  <>f  IFWC  for  Station  161770,  SKMI1W1  . \T1\SK 


121.  Cumulative  Frequency  Distribution  of  ILWC  for  Station  307580,  CHITA 
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Figure  A23.  Comparison  by  Station  of  1973  and  the  5-yr  Sample  in  Terms  of 
Average  ESI  for  Winter  and  Spring 
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Figure  A24.  Comparison  by  Station  of  1973  and  the  5-yr  Sample  in  Terms  of 
Average  ESI  for  Summer  and  Autumn 
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Figure  A25.  Comparison  by  Station  of  1973  and  the  5-yr  Sample  in  Terms  of 
Average  ILWC 


Table  A3.  Comparison  by  Station  and  Season  of  1973  and  the  5-yr  Sample  in  Terms  of  the  90th  Percentile  of  ESI.  (gkm 


Table  A5.  Comparison  by  Station  and  Season  of  1973  and  the  5-yr  Sample  in  Terms  of  Average  ILWC.  (g  m 


Comparison  by  Station  and  Season  of  1973  and  the  5-yr  Sample  in  Terms  of  the  90th  Percentil 


